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ABSTRACT 

■ We present XMM-Newton /EPIC spectra for the Laor et al. sample of Palomar Green 

^3 ' quasars. We find that a power-law provides a reasonable fit to the 2-5 keV region of 

t^. , the spectra. Excess soft X-ray emission below 2 keV is present for all objects, with the 

C ■ exception of those known to contain a warm absorber. A single power-law is, however, 

a poor fit to the 0.3-10.0 keV spectrum and instead we find that a simple model, 
consisting of a broken power-law (plus an iron line), provides a reasonable fit in most 
cases. The equivalent width of the emission line is constrained in just twelve objects 
but with low (< 2a) significance in most cases. For the sources whose spectra are well- 
ed . fit by the broken power-law model, we find that various optical and X-ray line and 
CL| continuum parameters are well-correlated; in particular, the power-law photon index 

is well-correlated with the FWHM of the H/3 line and the photon indices of the low 
£^ ■ and high energy components of the broken power-law are well-correlated with each 

^ | other. These results suggest that the 0.3-10 keV X-ray emission shares a common 

C$ . (presumably non-thermal) origin, as opposed to suggestions that the soft excess is 

directly produced by thermal disc emission or via an additional spectral component. 
We present XMM-Newton OM data which we combine with the X-ray spectra so as 
. to produce broad-band spectral energy distributions, free from uncertainties due to 

' long-term variability in non-simultaneous data. Fitting these optical-UV spectra with 

a Comptonized disc model indicates that the soft X-ray excess is independent of the 
accretion disc, confirming our interpretation of the tight correlation between the hard 
and soft X-ray spectra. 
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1 INTRODUCTION 

The launch of XMM-Newton has seen great advances in all 
areas of X-ray astronomy but its multiwavelength capabil- 
ities might be considered one of its greatest assets. This is 
of particular importance when it comes to studying quasars 
which are well-known for emitting at a broad range of fre- 
quencies. Different components of the quasar emit in differ- 
ent regions of the electromagnetic spectrum and so multi- 
wavelength analysis is vital. 

The "big blue bump" (BBB) which is observed in the 
optical-ultraviolet range, peaking at ~ 1000 A, is thought 
to be the signature of thermal emission from an accretion 
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disc (e.g. Peterson 1997 and references therein). Above 2 
keV the spectra can be well-modelled using a power-law, 
which is assumed to be due to Comptonization of lower en- 
ergy photons. The "soft excess" above the power-law, often 
seen in X-ray spectra of quasars at energies below ~ 2-3 
keV, has been attributed to the high energy end of the BBB 
(although see also Gierlihski & Done 2004); various models 
have been suggested for the origin of this soft X-ray excess 
such as direct thermal emission from the accretion disc or 
reprocessing of disc emission, (e.g. Peterson 1997 and ref- 
erences therein). Two additional components thought to be 
present in quasars, but not accessible to XMM-Newton, are 
the infrared bump (assumed to be emission from dust) and 
the radio jet (e.g. Peterson 1997 and references therein) . 

The Laor et al. (1994, 1997) complete sample of 23 
quasars was chosen from the Bright Quasar Survey (BQS), 



2 Brocksopp et al. 



Table 1. Table listing the 23 PG quasars in the Laor et al. sample, their redshifts and Galactic absorption and details of the XMM-Newton 
observations. Measurements of the Galactic absorption have been obtained from Laor et al. (1997 and references therein). 



Q uasar 


Rcdshift 




Observation 


XMM 


Observation 


Exposure 


Count Rate 






(xlO 20 cm" 2 ) 


Date 


Revolution 


ID 


(PN/MOS; ks) 


(PN, MOS1, MOS2; cts/s) 


PG 0947+396 


0.2060 


1.92 


2001-11-03 


349 


0111290101 


17.6/21.1 


1.68 0.31 0.31 


PG 0953+414 


0.2341 


1.12 


2001-11-22 


358 


0111290201 


10.6/13.2 


3.64 0.63 0.60 


PG 1001+054 


0.1610 


2.88 


2003-05-04 


623 


0150610101 


9.1/11.6 


0.04 0.01 0.01 


PG 1048+342 


0.1670 


1.74 


2002-05-13 


444 


0109080701 


28.1/32.2 


1.12 0.19 0.19 


PG 1114+445 


0.1440 


1.94 


2002-05-15 


445 


0109080801 


37.8/42.3 


0.66 0.18 0.17 


PG 1115+407 


0.1540 


1.74 


2002-05-17 


446 


0111290301 


15.0/20.3 


1.99 0.29 0.30 


PG 1116+215 


0.1765 


1.44 


2001-12-02 


363 


0111290401 


5.5/8.4 


4.79 0.74 0.73 


PG 1202+281 


0.1653 


1.72 


2002-05-30 


453 


0109080101 


12.9/17.4 


2.55 0.50 0.50 


PG 1216+069 


0.3313 


1.57 


2002-12-18 


554 


0111291101 


14.0/16.5 


1.03 0.19 0.18 


PG 1226+023 


0.1583 


1.68 


2000-06-13 


94 


0126700301 


44.3/58.5 


34.82 1.84 1.36 


PG 1309+355 


0.1840 


1.01 


2002-06-10 


458 


0109080201 


25.2/29.0 


0.43 0.08 0.08 


PG 1322+659 


0.1680 


2.88 


2002-05-11 


443 


0109080301 


8.7/11.6 


2.24 0.33 0.33 


PG 1352+183 


0.1520 


2.84 


2002-07-20 


478 


0109080401 


12.5/14.9 


2.19 0.34 0.34 


PG 1402+261 


0.1640 


1.42 


2002-01-27 


391 


0109081001 


9.1/12.0 


2.98 0.44 0.44 


PG 1411+442 


0.0896 


1.05 


2002-07-10 


473 


0103660101 


27.5/41.1 


0.11 0.02 0.02 


PG 1415+451 


0.1140 


0.96 


2002-12-08 


549 


0109080501 


21.2/24.1 


1.36 0.21 0.22 


PG 1425+267 


0.3660 


1.54 


2002-07-28 


482 


0111290601 


45.9/- 


0.63 0.45 0.45 


PG 1427+480 


0.2210 


1.69 


2002-05-31 


453 


0109080901 


35.2/39.0 


1.10 0.19 0.20 


PG 1440+356 


0.0790 


0.97 


2001-12-23 


373 


0107660201 


24.8/30.4 


6.26 0.78 0.77 


PG 1444+407 


0.2673 


1.09 


2002-08-11 


489 


0109080601 


18.7/21.5 


0.91 0.13 0.13 


PG 1512+370 


0.3707 


1.40 


2002-08-25 


496 


0111291001 


17.6/20.4 


1.45 0.28 0.29 


PG 1543+489 


0.4000 


1.68 


2003-02-08 


580 


0153220401 


8.2/9.8 


0.33 0.05 0.05 


PG 1626+554 


0.1330 


1.55 


2002-05-05 


440 


0109081101 


5.5/8.9 


2.92 0.44 0.47 



a subset of the Palomar Green (PG) survey compiled by 
Schmidt & Green (1983). The PG survey consists of 114 
AGN, 92 of which are quasars. With selection criteria Mb < 
—23 and U — B < —0.44, the BQS quasars are thought to 
be representative of most quasars (although see also Jester 
et al. 2005) ; since they were selected purely on their optical 
properties they were also thought to be unbiased in terms 
of their X-ray properties. In order to create a complete sam- 
ple, the Laor et al. sample was further constrained by two 
additional selection criteria: 

(i) z < 0.400, in order to prevent the rest frame 
0.2 keV being redshifted beyond the range observable by 
ROSAT/PSPC 

(ii) nfzcai < 1.9 x 10 20 cm~ 2 , in order to minimize the 
effects of Galactic absorption 

The full analysis of the ROSAT observations was pre- 
sented in Laor et al. (1997) and suggested that the 0.2-2 keV 
range of a quasar's rest frame spectrum could be fit by a sin- 
gle power-law, with no evidence for any soft excess emission 
and an upper limit (95% confidence) of ~ 5 x 10 19 cm~ 2 
on the amount of excess foreground absorption by cold gas. 
These data also suggested that the mean soft X-ray contin- 
uum agreed with an extrapolation of the mean UV contin- 
uum. The authors concluded that a thin bare accretion disc 
model was unable to reproduce the 0.2-2 keV spectral shape 
(Zheng et al. 1997; Laor et al. 1997). 

We have used XMM-Newton to update the ROSAT ob- 
servations, providing a much wider X-ray coverage (0.3-10 
keV), higher energy resolution and a high S/N ratio, and also 
providing simultaneous optical and UV data from the Op- 
tical Monitor (OM). Our aims were to investigate whether 



or not claims to the presence of a soft excess, intrinsic ab- 
sorption and/or iron lines were justified and physical when 
viewed in the context of the broad-band spectra. We have 
restricted our observations to the complete Laor et al. sam- 
ple of 23 quasars but compare our X-ray spectral parameters 
with those obtained for alternative BQS quasars by Porquet 
et al. (2004; their sample consisted of 21 PG quasars, 15 of 
which were Laor objects) and Piconcelli et al (2005) . In this 
paper we present the observations in Section 2, the X-ray 
spectroscopy in Section 3, line and continuum correlations 
in Section 4, the broad-band spectral energy distributions 
in Section 5 and discuss our results in Section 6. 



2 OBSERVATIONS 

XMM-Newton observed the 23 Laor targets between 2001 
November and 2003 May and obtained images with each of 
the three X-ray EPIC cameras (pn, MOS1, MOS2); details 
of the sources and their observations are listed in Table 1. 
Each of the three X-ray cameras has a spectral range of 
0.2-12 keV, spectral resolution of 20-50 and an angular res- 
olution of ~ 6". The EPIC large window mode and thin 
filter were used (with the exception of PG 1411+442 and 
PG 1425+267 which were observed in full frame mode, and 
PG 1226+023 which was observed with the medium filter 
and small window mode in order to minimise pile- up). The 
data were reduced using the SAS v5.4 epchain and em- 
CHAIN pipelines for the pn and MOS respectively. We re- 
reduced a selection of the sources with the SAS v6.0 version 
of the software but found it made no significant difference 
to the subsequent fits. All datasets were inspected for pile- 
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up by comparing the preliminary lightcurves with the pub- 
lished threshold count rates for each camera, as well as the 
epatplot routine; PG 1226+023 (= 3C 273) was the only 
source affected by pile-up. Response matrices and auxiliary 
files were generated using the rmfgen and ARFGEN tasks 
respectively. 

Longer wavelength images of each quasar were ob- 
tained, simultaneously with the X-ray data, using the Op- 
tical Monitor through the optical U, B and V filters and 
the UV UVW1, UVM2 and UVW2 filters (Mason et al. 
2001). They were reduced using the SAS v6.0 OMICHAIN 
reduction pipeline, which uses an aperture of 12 pixels. 
We list the central wavelength of the filters and resultant 
photometry in Table 2. Finally the magnitudes were cor- 
rected for Galactic extinction using the extinction law of 
Cardelli et al. (1989) and flux calibrated with respect to 
Vega. Flux conversions were obtained using the forumla 
Mvega - M quaaa r = -2.5 log(F V e g a / F qua sar ) , where Vega 
has a brightness of 0.03 magnitudes in the B and V fil- 
ters, 0.025 magnitudes in the U and UV filters and fluxes 
of 3.2652743, 5.9712835, 3.1594819, 3.7391873, 4.4609810, 
5.3884371 xl0~ 9 erg/cm 2 /s/A for the V, B, U, UVW1, 
UVM2 and UVW2 filters respectively 1 . 



3 X-RAY SPECTROSCOPY 

The task xmmselect was used to extract an X-ray spec- 
trum from a circular region centred on the source. Further 
circular regions were positioned on source-free areas of sky 
to provide a background sample. We corrected for the pile- 
up in PG 1226+023 by running the SAS v6.0 chains and 
extracting the spectrum from an annulus centered on the 
centre of the source. Photons with pattern 0-4 and 0-12, 
for the pn and MOS respectively, and with quality flag 
were included, grppha was used to bin the data using a 
minimum of 20 counts in each bin and 150 counts per bin in 
the case of PG 1626+554. The extracted spectra were anal- 
ysed using the xspec spectral fitting package (Arnaud et al. 
1996). 

Initial fits were made using a simple absorbed power- 
law, taking into account both Galactic absorption (obtained 
from Laor et al. 1997 and references therein) and possible 
redshifted absorption intrinsic to the quasar. We restricted 
the energy range to 0.3-10 keV in the observed frame for 
which the instrument response is generally well-known. It 
was clear however that there was a significant discrepancy 
between the pn and MOS spectra with as much as a 10% 
higher value of the power-law photon index for the pn rel- 
ative to that of the MOS. Kirsch (2003) outlines two pos- 
sible explanations for this problem; firstly a gradual change 
in the MOS response from revolution 200, and particularly 
450, onwards (as discussed by Porquet et al. 2004). The sec- 
ond factor is a ~ 10% higher flux for the pn in the 0.3- 
1.0 keV range compared with the MOS, which is possibly 
due to uncertainties in the vignetting and CCD-quantum- 
efficiency. Since our sources are all relatively soft we suggest 
that both factors are contributing to the discrepancy. While 

l 
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the first factor is an uncertainty with the MOS calibration, 
we note that the second factor is not necessarily restricted 
to the MOS. We therefore do not discard all the MOS data 
(contrary to Porquet et al. 2004) but discard only the MOS 
data below 0.7 keV. We find that the pn and MOS are in 
good agreement (typically to within 90% confidence con- 
tours) above this energy. We re-reduced four sources chosen 
at random using the SAS v6.0 software and found that there 
was no significant improvement in pn/MOS agreement. 

We then performed power-law fits to the spectra, fit- 
ting pn and MOS simultaneously and restricting the energy 
range to 2-5 keV in the rest frame (in contrast to Porquet 
et al. 2004). We included only the effect of Galactic absorp- 
tion and note that this region of the spectrum is useful for 
preliminary fits, since it is typically thought to be free of 
model components due to soft excess, absorption, reflection 
or emission lines. A power-law was indeed a good fit to the 
spectra in this range and the resultant parameters are listed 
in Table 3. The values of F fall in the range 1.4-2.6, with 
the exception of PG 1001+054 and PG 1411+442. Extrap- 
olating the 2-5 keV fit to the full XMM range confirms the 
result of Porquet et al. (2004) that a soft excess is present 
in all cases, except for those complicated by the presence of 
a warm absorber (Ashton et al. 2004; Ashton, Brocksopp et 
al. in prep.). 

Indeed, over the full XMM-Newton energy range (0.3- 
10 keV) a simple absorbed power law is not a good fit to 
the data (see xi m final column of Table 3). Many of the 
sources show excess emission in the low and/or high en- 
ergy regions of the spectrum. Again, this confirms the re- 
sults of Porquet et al. (2004) and Piconcelli et al. (2004) 
and so we do not repeat the work here. Instead we try fit- 
ting a broken power-law, again including Galactic absorp- 
tion and possible intrinsic, redshifted absorption (Fig. 1). 
This gives a more acceptable fit (0.94 < xl < 1-13), with 
the exception of PG 1114+445, PG 1216+069, PG 1226+023 
and PG 1411+442 as discussed below (Table 4). Intrinsic 
broad-band absorption is not required in the majority of 
our sources, although Porquet et al. (2004) included narrow 
absorption edges for some objects. Values of the photon in- 
dex for the low energy power-law fall in the 2.0-4.0 range; 
for the high energy power-law the photon indices lie in the 
range 1.5-2.9 (again with the exceptions of PG 1001+054 
and PG 1411+442). The break energy was in the range 0.7- 
2.1 keV (or 0.9-2.5 keV in the rest frame). 

We also include a single iron line of width 0.01 or 1.0 
keV in our model, placing it at each of the three energies 
6.4, 6.5 and 6.7 keV (Table 5). However, we find that the 
equivalent widths of any potential narrow or broad lines can 
be constrained in just twelve sources; with most of these de- 
tections having a significance of < 2a, it is debateable as to 
whether these lines can be considered real. 99% confidence 
upper limits to the equivalent width are also listed in Ta- 
ble 5. The central energies of the detected lines are typically 
consistent with 6.4 or 6.5 keV, corresponding to neutral and 
weakly ionized iron emission. A more detailed search for iron 
emission in PG 1226+023 has been published by Page et al. 
(2004a) and a weak broad line found by co-adding spectra. 
We note that Jimenez-Bailon et al. (2005) study a sample 
of 38 PG quasars and detect iron lines in 20 of them, 13 of 
which are included in the Laor et al. sample. This does not 
necessarily contradict our results as these authors use a va- 
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Table 2. Table of optical and UV photometry from the Optical Monitor observations. The second column under each filter heading gives 
the statistical error. Details of the OM and the filters are given in Mason ct al. (200f). Flux conversions are stated in the text. 



Filter/magnitudes (A e ff) 

Quasar V (5483 A) B (4443 A) U (3735 A) UVW1 (2910 A) UVM2 (2310 A) UVW2 (2120 A) 



PG 


0947+396 










15.31 


0.01 


15.13 0.01 


14.92 


0.01 


14.77 0.01 


PG 


0953+414 














13.16 0.01 


13.71 


0.01 


13.55 0.01 


PG 


1001+054 






















PG 


1048+342 


16.51 


0.03 


16.80 


0.02 


15.66 


0.01 


15.44 0.01 


15.36 


0.02 


15.25 0.02 


PG 


1114+445 


15.78 


0.01 


16.17 


0.01 


15.08 


0.01 


15.00 0.01 


15.15 


0.01 


15.08 0.02 


PG 


1115+407 






16.00 


0.01 


14.83 


0.01 


14.54 0.01 


14.25 


0.01 


14.12 0.01 


PG 


1116+215 


- 




14.32 


0.01 






12.85 0.01 






12.43 0.01 


PG 


1202+281 


16.27 


0.03 


16.52 


0.02 


15.32 


0.01 


15.10 0.02 


14.98 


0.02 


14.80 0.03 


PG 


1216+069 


15.43 


0.01 


15.64 


0.01 


14.60 


0.01 


14.46 0.01 


14.20 


0.01 


14.13 0.01 


PG 


1226+023 


12.64 


0.01 


12.94 


0.01 


11.77 


0.01 


11.49 0.01 


11.33 


0.01 


11.28 0.01 


PG 


1309+355 


15.64 


0.01 


16.03 


0.01 


15.00 


0.01 


14.84 0.01 


14.76 


0.01 


14.70 0.02 


PG 


1322+659 










14.82 


0.01 


14.58 0.01 


14.36 


0.01 


14.27 0.02 


PG 


1352+183 


16.37 


0.03 


16.63 


0.02 


15.41 


0.01 


15.14 0.01 


14.98 


0.02 


14.85 0.03 


PG 


1402+261 














14.04 0.01 


13.74 


0.01 


13.61 0.01 


PG 


1411+442 










13.69 


0.01 










PG 


1415+451 


15.96 


0.02 


16.54 


0.01 


15.43 


0.01 


15.18 0.01 


15.15 


0.02 


15.02 0.02 


PG 


1425+267 


17.16 


0.04 


17.17 


0.02 


16.05 


0.01 


16.10 0.01 


15.90 


0.02 


15.79 0.03 


PG 


1427+480 


16.63 


0.04 






15.69 


0.01 


15.53 0.01 


15.40 


0.01 


15.22 0.02 


PG 


1440+356 


14.63 


0.01 






13.92 


0.01 


13.66 0.01 


13.54 


0.01 


13.43 0.01 


PG 


1444+407 


15.94 


0.02 


16.07 


0.01 


14.95 


0.01 


14.81 0.01 


14.48 


0.01 


14.39 0.02 


PG 


1512+370 


16.65 


0.03 


16.75 


0.01 


15.67 


0.01 


15.58 0.01 


15.06 


0.02 


14.85 0.02 


PG 


1543+489 










16.50 


0.01 




15.14 


0.01 




PG 


1626+554 


15.79 


0.02 










14.60 0.01 


14.49 


0.01 


14.43 0.02 



Table 3. Table of parameters resulting from a single power-law fit to the 2—5 keV energy range (rest frame; column 2 lists the equivalent 
observed frame energies), including only Galactic absorption. The penultimate column lists the luminosity of each quasar, assuming a 
standard WMAP cosmology. Clearly a power law spectrum is a reasonable assumption in this region of the spectrum for the majority of 
the sources. The final column shows the reduced x 2 f° r a power-law model fit to the full 0.3-10 keV range. Clearly this is a poor fit in 
all cases. 



Quasar 


Observed Frame Photon Index 


X 2 (DOF) 


x 2 


Probability 


Luminosity 


xl 




Energy (keV) 








(xlO 44 ergs" 1 ) 


(0.3-10 keV) 



PG 0947+396 


1.7 - 


4.1 


2 10+ ' 08 

z - ±u -0.08 


279.89 


(286) 


0.98 


0.59 


1.18 


1.38 


PG 0953+414 


1.6 - 


4.1 


2 19+ -" 7 
z - ±M -0.07 


302.43 


(290) 


1.04 


0.30 


2.89 


1.65 


PG 1001+054 


1.7 - 


4.3 


14+ - 46 


1.14 


(8) 


0.14 


1.00 


0.025 


4.02 


PG 1048+342 


1.7- 


4.3 


1 96+ 09 

^ -0.09 


194.69 


(213) 


0.91 


0.81 


0.55 


1.45 


PG 1114+445 


1.7 - 


4.4 


1 36+ 006 
1 - JD -0.05 


412.14 


(443) 


0.93 


0.85 


0.57 


4.64 


PG 1115+407 


1.7 - 


4.3 


2 o 4 +0.l0 
z -° -0.10 


170.59 


(192) 


0.89 


0.86 


0.50 


1.55 


PG 1116+215 


1.7 - 


4.2 


9 Q9+0.09 
z - zz -0.09 


226.40 


(217) 


1.04 


0.32 


1.79 


1.51 


PG 1202+281 


1.7 - 


4.3 


1 76+ ' 07 
1 -'°-0.07 


290.70 


(342) 


0.85 


0.98 


1.36 


1.54 


PG 1216+069 


1.5 - 


3.8 


2 02+ ' 12 
z - uz -0.12 


184.09 


(142) 


1.30 


0.01 


2.14 


1.97 


PG 1226+023 


1.7 - 


4.3 




1425.47 


(864) 


1.65 


le-13 


27.0 


4.99 


PG 1309+355 


1.7- 


4.2 


l.80tS;ll 


121.05 


(137) 


0.88 


0.83 


0.34 


1.74 


PG 1322+659 


1.7- 


4.3 


2 o 4 +o.i3 

z -° -0.13 


120.17 


(117) 


1.03 


0.40 


0.64 


1.48 


PG 1352+183 


1.7 - 


4.3 


2 10+ ' 11 


142.59 


(146) 


0.98 


0.56 


0.66 


1.39 


PG 1402+261 


1.7 - 


4.3 


2 42+ - 10 
z -* z -o.io 


185.92 


(189) 


0.98 


0.55 


0.82 


1.63 


PG 1411+442 


1.8 - 


4.6 


n 34+019 
U - J4 -0.31 


46.84 


(40) 


1.17 


0.21 


0.02 


7.33 


PG 1415+451 


1.8 - 


4.5 


2 14 +0.09 
z - ± ^-0.09 


174.80 


(195) 


0.90 


0.85 


0.21 


1.68 


PG 1425+267 


1.5 - 


3.7 


1 -^ a -0.07 


436.24 


(379) 


1.15 


0.02 


3.05 


1.47 


PG 1427+480 


1.6 - 


4.1 


2 04+ ' 07 
z - u4 -0.07 


296.92 


(327) 


0.91 


0.88 


0.89 


1.45 
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Figure 1. A plot for each quasar showing the X-ray spectra (pn and MOS) in ^F„-space and the best fit to a model comprising an 
absorbed (allowing for both Galactic + intrinsic absorption) broken power-law plus iron line. This simple model provides a good fit to 
the spectra for the majority of the sources. We note that Ashton et al. (2004) have also included a warm absorber component in order to 
fit PG 1114+445 and PG 1309+355 successfully and suspect that similar detailed modelling would be required in order to fit the spectra 
of PG 1425+267 and PG 1411+442 (Ashton, Brocksopp et al. in prep.). 
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Table 4. Model-fit to the 0.3-10 keV (observed frame) data using an absorbed broken power-law with the addition of a single iron line 
of width 0.01 keV in the 5.5-8.0 keV (rest frame) range. The observed luminosity in this energy range has also been listed. This model 
provides a reasonably good fit for the majority of sources with the most notable exceptions being PG 1114+445 (see Ashton et al. 2004) 
and PG 1411+442 (which is unconstrained; sec Ashton, Brocksopp et al. in prep). 



Quasar 


Pholnd. 1 


Pholnd. 2 


Break E 


Absorption 


Cont. Norm. 


Luminosity 


X 2 (DOF) 


xl 


Probability 










(keV) 


(xlO 20 cm" 2 ) 


(xl0~ 3 ) 


(xlO 44 ergs s- 1 ) 










0947+396 


9 4 q+0.07 
^•^°-0.03 


1 07+0.13 
1 - 8 '-0.09 


9 no +0. 30 
2.02_ 044 


< 0.63 


0.84 


5.4 


786.29 (790) 


1.00 


0.53 


Dp 


uyoo+414 


z - uo -0.03 


9 n7 +0.06 
z - u '-0.08 


1 Kc+0.22 
L56 -0.19 


< 1.54 


1.69 


15.7 


826.15 (777) 


1.06 


0.11 


pp 

I V_J 


1UU1TUD4 


A Qfi + 172 
4 - db -0.54 


n 74 +0.40 
U -' 4 -0.39 


1 n o+0.21 


0.00 


n m 

U.Ul 


u. 1 


97 7Q 




U.O / 


pn 
r v_i 


1 C\A Q 1 Q A O 

!U4o+o4z 


9 OQ + 0.11 


1 - 84 -0.09 


1 fifi+ ' 30 
1-DO_0.14 


1 cj; 


U.OD 


n 


71 A 1 K f"£;7c^ 

( 14. ID 1^0 ( J 


l.UD 


n 1 a 

U. 14 


pp 


1 ~l 1 A J_A AK 
1 1 14+440 




















pp 


111 Z,J_AC\7 


O QQ+ - 17 

z - ya -o.n 


9 o fi +0.06 


, oc+0.22 
J^-^O-O.IS 


z. / 


u.yu 


0.0 


uucii ^U4:y^ 




U.O / 


PG 


1116+215 


z -'O_0.04 


9 iq+0.08 
z ' la -0.08 


1 4R+0-22 
i - 4t >_0.22 




9 it 
z. 10 


1 


709 S 1 ^ ^77^ 


1 04 


24 


PG 


1202+281 


9 97 +0.04 
z - z '-0.03 


1 71+O.O5 
1 -' L -0.05 


1 70+°' 14 

1- < u -0.17 


< 0.50 


1.24 


5.5 




0.98 


0.61 


PG 


1216+069 


3 r)9+ ' 32 


1 O7+0.08 
-0.08 


96 +0 17 


4 25+ 7 - 55 


0.39 


2.3 


655.40 (480) 


1.37 


2c-7 


PG 


1226+023 


2.04^1 


l-68±g : gl 




< 0.12 


21.52 


94.4 


4353.74 (2785) 


1.56 


<lc-13 


PG 


1309+355 




1 oq+0.07 
1 - o,:> -0.06 




6.77t 4 ;™ 


0.13 


1.3 


456.39 (473) 


0.96 


0.70 


PG 


1322+659 


O QO+0.19 

°- uo -0.09 


9 91+O.O9 
z ' z -0.12 


1 49+ 015 


< 2.76 


1.00 


4.3 


476.39 (495) 


0.96 


0.72 


PG 


1352+183 


2 61+ 05 


1 99+ 011 
^ -0.07 


1.91+0' 24 

— 0.25 


< 0.55 


1.07 


3.5 


587.43 (561) 


1.05 


0.21 


PG 


1402+261 


9 07+O.II 
z - 8 '-0.04 


9 90+O.O9 

^•^-o.io 


1 cc+0.23 
± - oo -0.14 


< 1.31 


1.31 


5.4 


727.0 (642) 


1.13 


0.01 


PG 


1411+442 


07+O.47 

•^'-o.is 


S5+ 011 


1 34+ 009 


< 0.04 


0.02 


0.1 


308.3 (157) 


1.96 


6e-12 


PG 


1415+451 


9 ro+0.08 
z ' oa -0.03 


2 03+ 09 
z ' uo -0.07 


1 66+ 016 

1 - oo -0.21 


< 0.79 


0.62 


1.2 


698.40 (671) 


1.04 


0.23 


PG 


1425+267 


oq+0 93 
°- ua -0.39 


1 4Q+0.04 


67+ 05 


4 fiQ+ 8 ' 33 
^• D9 -3.31 


0.16 


9.9 


438.37 (409) 


1.07 


0.56 


PG 


1427+480 


2 48+ ' 04 
z '^°-0.03 


1 96+ 07 
J -- au -0.07 


1 oi+0.20 
1 - fil -0.26 


< 1.27 


0.55 


4.2 


821.74 (843) 


0.97 


0.69 


PG 


1440+356 


10+O.12 
°- lo -0.09 


9 40+0.07 


1 43+0.13 


1 - iO -0.79 


2.50 


2.2 


1007.30 (921) 


1.09 


0.02 


PG 


1444+407 


2 84 +0 11 


9 Q7+0.21 
z ' u '-0.20 


2 00+ ' 37 
z ' uu -0.46 


< 1.34 


0.40 


5.2 


480.77 (460) 


1.05 


0.24 


PG 


1512+370 


9 To+0.18 
z -°°-0.04 


1 otr+0.09 
l-»O_0.06 


1 41+0-21 
± -^ ± -0.34 


< 1.30 


0.71 


19.1 


679.47 (663) 


1.02 


0.32 


PG 


1543+489 


n 07+1.79 


9 ocr+0.34 


90+ 2 ' 25 
u ' au -0.20 


< 32 


0.15 


4.6 


164.93 (158) 


1.04 


0.34 


PG 


1626+554 


2 44+ - 24 


9 nt -+0.06 
z - uo -0.13 


1 n rr+0.57 
1 - uo -0.16 


< 2.24 


1.36 


3.6 


180.47 (166) 


1.09 


0.21 



riety of different models, some of which we show in Sections 
4-6 may be inappropriate. Further discussion of the iron 
lines in these sources, with reference to reflection models, 
can also be found in Porquet et al. (2004). 
Alternative 

models, such as black-body, bremsstrahlung, Comptoniza- 
tion and more complicated absorption components, are of- 
ten used to improve fits. Porquet et al. (2004) and Piconcelli 
et al. (2004) have performed detailed fitting of these mod- 
els to sub-sets of the Laor et al. sample and so we do not 
repeat their work. These authors found that the addition of 
black-body or bremsstrahlung model components improved 
fits but that no single model was able to fit all the objects. 
The values of the resultant parameters were also regarded 
as phenomenological rather than physical, in particular re- 
quiring temperatures much higher than those thought to be 
reached in an accretion disc (Piconcelli et al. 2004). Comp- 
tonization models tended to be more successful and provided 
a good fit to the soft excess of the PG quasars (Porquet et al. 
2004; Piconcelli et al. 2004; Gierlinski & Done 2004). While 
the Comptonization of disc photons is indeed a potential 
physical model, there are still problems in that the resultant 
temperatures lie in an unrealistically narrow range (0.1-0.2 
keV), which does not reflect the wider range of probable 
disc temperatures, and the soft X-ray emitting region would 
have to be extremely compact (R ~ 10 12 cm). 

The few sources which are not well-described by 
the simple broken power-law model merit further discus- 
sion. Ashton et al. (2004) perform detailed modelling of 



PG 1114+445 and PG 1309+355 and require the addition of 
a warm absorber component. PG 1425+267 is qualitatively 
similar to PG 1309+355 and may also require an additional 
absorption component (as suggested previously by Reeves 
& Turner 2000). PG 1001+054 and PG 1411+442 are sig- 
nificantly different from the other sources; PG 1411+442 
in particular appears to require an additional high energy 
component, perhaps reflection which would be consistent 
with the apparent strong emission line present in this source. 
Broken power-law fits to the spectra of PG 1216+069 and 
PG 1226+023 are also unacceptable, although in these cases 
it is not so clear why; the PG 1226+023 spectra are of much 
higher S/N than the other sources and this seems to be the 
reason why the fit is not as successful (suggesting, perhaps, 
that higher S/N observations of the other sources would also 
require more sophisticated models). Detailed modelling of 
individual sources is beyond the scope of this paper but will 
be addressed in follow-up work. 



4 X-RAY/OPTICAL CONTINUUM AND LINE 
LUMINOSITY CORRELATIONS 

We have plotted various combinations of spectral parame- 
ters and line and continuum 2 luminosities in Figures 2 and 
3. The optical data were obtained from Laor et al. (1997) 

2 Here we use only the luminosity derived using the pn data since 
the MOS spectra were unreliable at 0.3 keV 
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Figure 2. Scatter plots showing various optical parameters (obtained from Laor et al. (1997)) plotted against the broken power-law 
parameters obtained in Section 3 (omitting PG 1114+445). The break energies show no indication of correlated behaviours; however, 
with the exception of PG 1001+054 and PG 1411+442, the [OIII] luminosity and, particularly, the H/3 FWHM are correlated with the 
photon index, in the sense that as the spectrum hardens the [OIII] luminosity and H/3 FWHM increase. Spearman rank order correlation 
coefficients for each plot are listed in Table 6. 
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Table 5. Equivalent widths of a single possible iron line, of width 0.01 or 1.0 keV and rest frame energies 6.4, 6.5 or 6.7 keV, of an 
absorbed broken power-law model. Upper limits on the equivalent width are to within 99%. We note that Page ct al. (2004a) detected a 
weak, broad line in PG 1226+023 on co-adding a number of different observations. ^The equivalent width for PG 1114+445 is obtained 
from Ashton et al. (2004) who include a warm absorber in the model. 
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but we note that the lack of simultaneous observations of 
any varying sources is more likely to weaken the correla- 
tions than to strengthen them. We then calculated the as- 
sociated Spearman rank order correlation coefficients (p); 
many of the plots show only scatter, but there are also a 
number of significant correlations (0.6 < p < 0.9 to within 
99% confidence). We list these higher values of p in Ta- 
ble 6. We note that PG 1114+445 has been omitted from 
the plots on account of its extremely poor fits to the bro- 
ken power-law model. For each pair of parameters in the 
table, the first value of p is for the plots from Figures 2 and 
3; the second value is obtained through further omission of 
PG 1001+054, PG 1216+069, PG 1411+442, PG 1309+355 
and PG 1425+267, sources which have been labelled "X-ray 
weak" by Laor et al. (1997), shown to contain a warm ab- 
sorber (e.g. Ashton et al. 2004; Reeves & Turner 2000) or 
for which the broken power law provides an inadequate fit 
to the spectrum. 

Allowing for the few discrepant sources which require a 
more complicated model, it is clear that there are a num- 
ber of strong correlations present. In Figure 2 we show that 
both the hard and soft photon indices are correlated with 
the [om] luminosity and the H/3 FWHM . The break en- 
ergies do not show any form of correlated behaviour with 
other quantities. Similarly the 3000 A and H/3 luminosities 
do not seem to show any dependency on the X-ray spec- 
tral parameters. Figure 3 shows that the 2.0 and 10.0 keV 
luminosities are clearly correlated with all four optical pa- 
rameters; the 0.3 keV luminosities are not correlated with 
any of the optical parameters. 



These results confirm the results of Laor et al. (1997), 
Reeves & Turner (2000) and Brandt, Mathur & Elvis (1997) 
who find similar correlations, most notably the result that 
the X-ray spectrum steepens as the FWHM of the H/3 
line decreases. Reeves & Turner (2000) suggest that, since 
objects with narrower H/3 lines are thought to accrete at 
a higher fraction of the Eddington limit, there is greater 
Compton cooling of the hard X-ray-emitting corona and 
thus a steepening of the spectrum. They also suggest that 
a soft excess component associated with the accretion disc 
should therefore be present in the objects with narrower 
lines. However, we find that the soft excess is actually 
present below 2 keV in all sources in the Laor et al. sam- 
ple which do not contain some sort of warm absorber. In- 
stead we compare the photon indices of the low and high 
energy power-laws used in our broken power-law fits and 
find a strong correlation, confirming the result of Porquet 
et al. (2004); p = 0.97 (to within > 99% confidence) is ob- 
tained when we omit the (consistently discrepant) sources 
PG 1001+054, PG 1411+442, PG 1216+069, PG 1309+355 
and PG 1425+267. We note that the correlation is stronger 
for our sample than for that of Porquet et al. (2004), who ob- 
tain p = 0.85; this may be a reflection either on the slightly 
different sample or on our use of the broken power-law, in- 
stead of fixed energy ranges in the observed frame. We dis- 
cuss the implications of this correlation further in Section 
6. 
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Figure 3. Scatter plots showing various optical parameters (obtained from Laor et al. (1997)) plotted against X-ray luminosities at 
0.3, 2.0 and 10.0 keV (omitting PG 1114+445). It is clear that as the 2.0 and 10.0 keV luminosities increase, so too do the values of 
all four optical parameters; any correlation at 0.3 keV is less convincing, perhaps on account of the EPIC response. Pearson correlation 
coefficients for each plot are listed in Table 6. 



5 SPECTRAL ENERGY DISTRIBUTIONS 

The most comprehensive compilations of AGN / quasar spec- 
tral energy distributions (SEDs) presented to date are those 
of Elvis et al. (1994) and Kuraszkiewicz et al. (2003). These 



authors included 2-10 keV X-ray data obtained by the 
HEAO missions, archive data obtained from the literature 
and Kuraszkiewicz et al. (2003) also presented far-infrared 
ISO data. In particular, the Kuraszkiewicz et al. (2003) sam- 
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Table 6. Spearman rank order correlation coefficients (p > 0.5 to 99% confidence) calculated for each of the scatter plots in Figures 2 
and 3. Optical/UV luminosities arc from Laor ct al. (1997) and converted into the WMAP equivalent. For each pair of parameters, 
the first value of p is for the plots from Figures 2 and 3. The second value is obtained through further omission of the five discrepant 
objects; these quasars may be discrepant on account of being an "X-ray weak quasar" and containing a warm absorber (PG 1001+054, 
PG 1309+355, PG 1411+442 and PG 1425+267.) or on account of being less-well fit by a broken power-law than the majority of other 
sources (PG 1216+069). We note that a number of these correlation coefficients can be increased still further through omission of only 
PG 1001+054 and PG 1411+442, but do not include these values for consistency with Fig. ??. 



Luminosity (3000 A) 


Luminosity (H/3) 


Luminosity ([oiii]) 


FWHM (H/3) 


Low Energy T 






-/-0.66 


-0.53/-0.86 


High Energy T 






-0.52/-0.60 


-0.52/-0.79 


Break Energy 










Luminosity (0.3 keV) 


0.59/0.69 


0.52/0.48 






Luminosity (2.0 keV) 


0.77/0.77 


0.82/0.73 


0.63/0.65 


0.55/0.45 


Luminosity (10.0 keV) 


0.67/0.58 


0.78/0.66 


0.74/0.78 


0.78/0.73 


Low Energy T vs. High Energy r: p ~ 0.97 to 99% confidence (see text and Fig. 4) 



In* 



PG 1216+069 



PG 1425+267 



PC '411 +442 I 



PG 1543+489 



PG 1309+355 



PG 1001+054 



3 4 
Low Energy Photon Index 




2.5 3.0 
Low Energy Photon Index 



3.5 



Figure 4. LHS: Scatter plot showing the high energy photon indices plotted against the low energy photon indices. Clearly, with 
the exception of the six discrepant objects, there is a strong correlation (p ~ 0.97 to 99% confidence). This suggests that the regions 
responsible for the emission are either one and the same or in direct "contact" with each other. The discrepant objects are mainly those 
which contain warm absorbers (Ashton et al. 2004) and/or arc defined as "X-ray weak" (Laor et al. 1997). PG 1216+069 is poorly-fit by 
the broken power-law but it is less clear why it should be discrepant. RHS: Expanded view of the highly-correlated region of the plot. 
The dotted lines correspond to Tjjig^ = r^ ow and Vjjig^ = 0.39 + 0.63r£ otl ,, where the errors on the intercept and gradient are 0.08 
and 0.04 respectively. 



pie is hard X-ray selected which is unbiased by effects due 
to obscuration and reprocessing. 

The Laor et al. (1997) sample was optically-selected 
with the intent of defining a representative sample of X- 
ray properties; however we note that this has the effect of 
biasing the sample against obscured sources, as discussed 
by Kuraszkiewicz et al. (2003). Nonetheless, the quality of 
the XMM-Newton data that we present here is sufficiently 
superior to the HEAO data that compiling SEDs is a valu- 
able exercise, even if the multiwavelength results cannot be 
considered representative of all quasars. 

In Fig. 5 we plot the optical-X-ray spectra (OM pho- 
tometry plus EPIC spectra) for each quasar in the Laor et 
al. sample. They have been plotted in vF v space in Fig. 5 
and have been corrected for Galactic absorption. Clearly the 



plots show a wide range of spectral shapes and we discuss 
this in detail below. 

The SEDs are dominated, as expected (e.g. Elvis et al. 
1994), by the BBB in the optical-UV region. The turnover 
of the BBB is not observed for any of the Laor et al. sam- 
ple. However, as already noted, the majority of the X-ray 
spectra show a soft excess. Analysis of ROSAT data by 
Puchnarewicz et al. (1996 and references therein) revealed 
evidence for correlations between the optical and X-ray spec- 
tral slopes and the optical:X-ray flux ratio. Such correlations 
led to the conclusion that there is a strong behavioural link 
between the BBB and the soft X-ray excess. The additional 
correlations reported in Section 4 and references therein sup- 
port this conclusion. The correlation between the soft and 
hard X-ray power-law slopes suggests further still that there 
must be a behavioural link between the BBB and the X-ray 
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spectrum as a whole, as opposed to the soft excess in isola- 
tion. 

We have fit our XMM-OM photometry with Comp- 
tonized accretion disc curves generated by a modified Czerny 
& Elvis (1987) model (Siemiginowska et al. 1995; Starling 
et al. 2004). The accretion disc spectra are calculated as- 
suming that the disc is geometrically thin, emitting locally 
as a modified blackbody (including electron scattering and 
Comptonization in the atmosphere of the disc) in a Kerr ge- 
ometry. They do not include the effects of irradiation but are 
appropriate for simple modelling of the soft excess of X-ray 
spectra. The best-fit models are plotted on top of the data 
in Fig. 5 (solid lines) and appear to be a good or reasonable 
representation of the OM data for a number of objects. The 
values of xt are listed with the best-fit masses in Table 7; 
while some values of xt are relatively high we note that 
for this preliminary study of the SEDs we have used only 
a conservative estimate (optical: 1%, UV: 5%) for the flux 
calibration systematic errors and have relatively few points 
available for constraining the fits. The results also include 
estimates of the mass accretion rate and inclination angle; 
the full range of parameter-space is sampled in each case and 
the results also listed in Table 7. We inclulde approximate 
error ranges (to within 99% confidence) for the masses but 
were unable to constrain the mass accreion rate, luminos- 
ity or inclintaion angle. In order to determine the errors we 
have obtained the volume in parameter-space that is con- 
tained within Ax of the minimum \ for each fit. These 
error ranges should, however, be used with caution since 
they do not reflect exactly those points of parameter-space 
which were included within the A% 2 + Xmm volume. 

It is clear from Fig. 5 that, despite reasonable fits to 
the OM data, the Kerr models often predict an excess of 
flux compared with the observed X-ray data (6 out of 22 
cases). While we might expect this for those sources known 
to contain warm absorbers, this is not a satisfactory expla- 
nation for the remainder (e.g. PG 1115+407). If instead we 
use the published masses from the literature (obtained by 
line width techniques; Shields et al. 2003; Vestergaard 2002; 
Baskin & Laor 2005), the fits are much worse (>> 3a in 
most cases) and the predicted soft X-ray flux exceeds that 
observed in 7 systems. Most of the masses we derive, by 
fitting the OM data, are higher than published estimates 
in the literature; accretion rates in excess of the Eddington 
limit (and outside the model limits) and/or the effects of 
e.g. irradiation may be required in order to provide accept- 
able fits to these published mass values. In 11 objects (e.g. 
PG 0953+414), the models that use the published masses 
(dotted lines in Fig. 5) underestimate the optical/UV flux 
compared with the OM points. This is somewhat worrying 
since it shows that either the models are inadequate and e.g. 
irradiation needs to be included, or it is incorrect to assume 
that there is no significant contribution to the optical flux 
from e.g. host galaxy contamination and/or emission lines. 
Alternatively it may suggest that the source luminosity has 
varied since the mass estimates were published, as would be 
consistent with the discrepancies between the OM and IUE 
data in some objects. Finally, we do not confirm the corre- 
lation between mass and 2-10 keV luminosity (p ~ 0.5 to 
within 99.7% confidence) as found by Piconcelli et al. (2004), 
using either our derived masses or the published masses of 
the Laor PG quasars. 



The most striking feature of these plots is that there 
appears to be no relationship between the fits and the X- 
ray spectra. There is no source for which the soft excess 
has a sufficiently steep spectrum that it can be reproduced 
by the disc model, in contrast with a source like Mkn 841 
for which a disc spectrum can be fit to the combined soft 
excess and UV spectrum (Arnaud et al. 1985). Indeed, we 
find that the slope and flux of the soft excess is not predicted 
by the disc model to any significant degree, which provides 
an interesting challenge for interpreting the optical/UV/X- 
ray correlations discussed above. 

As an alternative, we have also performed fits to the 
spectra using Schwarzschild models, both with and without 
modifications to correct for the effects of opacity. The values 
of xt marked with "m" in Table 7 indicate which of the best- 
fit Schwarzschild models had been modified to include the 
effects of opacity. These fits are shown in Fig. 5 (dashed 
lines). Schwarzschild models do not predict significant flux 
in the X-ray region, and thus the problem of overestimating 
the X-ray flux is eliminated. Table 7 shows that the fits to 
the optical/UV data are improved for the majority of the 
Laor objects, compared with the Kerr fits (Ax 2 > LOcr for 
17 objects and A^ 2 > 2.6<r for 11 objects). There is also a 
marginal increase in agreement between the masses derived 
and the published values (dot-dashed lines in Fig. 5). 

Similar modeling of sub-sets of the Laor et al. sample 
was performed by other authors; we find that our resultant 
masses are comparable with, although better-constrained 
than, those of Laor (1990) but are consistently lower than 
those of Sun & Malkan (1989). We also note that Laor (1990) 
found that Kerr models were more successful at describing 
the data. Whilst there is evidence for a rotating black hole 
in MCG-6-30-15 (Wilms et al. 2001; Vaughan & Fabian 
2004) and broad iron lines are seen in several more AGN 
for which a relativistic disc line is a viable explanation (e.g. 
Mkn 766, Mason et al. 2003; Q0056-363, Porquet & Reeves 
2003; Mkn 205, Reeves et al. 2001), our preliminary study 
of the broad-band SEDs of Laor PG quasars might suggest 
that a Schwarzschild geometry is more appropriate. In any 
case the soft X-ray excess appears to be independent of the 
disc model. 



6 DISCUSSION 

The BQS is a well-studied sample and so there are various 
results in the literature with which to compare the XMM- 
Newton spectra. In particular, Laor et al. (1997) obtained 
values of the photon index for the 1.2-3.0 keV range which 
are comparable with the values we obtain for the softer of the 
two power laws, i.e. V — 2.2-3.7 below the break energy. We 
note that our derived break energies are in the range 0.9-2.6 
keV, whereas in effect the Laor et al. study was equivalent to 
assuming a break energy >2 keV, and so this a reasonable 
but not exact comparison 3 . Comparison of the XMM and 
ROSAT photon indices is shown in the top panel of Fig. 6. 

3 For simplicity and comparison with the literature, we discuss 
the spectra in terms of their power-law parameters but note that 
the successful fitting of e.g. blackbody and Comptonization mod- 
els to the soft X-rays by other authors may indicate a curved 
spectrum instead of a power-law. 
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Table 7. Results from fitting the Kerr and Schwarzschild SED models to the OM photometry. The first two columns list the sources 
and the number of degrees of freedom for the fits. The next four columns (first line of data) show the best-fit mass and x 2 of the 
Kerr and Schwarzschild models respectively. The seventh column lists the published mass-estimates for each source (from "Shields et al. 
2003, & Vcstcrgaard 2002; c Baskin & Laor 2005) and the final two columns show the minimum x 2 obtainable for these published masses. 
Derived masses for which 99% confidence intervals include the published mass values are marked by an asterisk. The values of x 2 marked 
with "m" indicate which of the best-fit Schwarzschild models had been modified to include the effects of opacity. We find that the values 
of xi are relatively high in some cases, but the majority appear to provide reasonably successful representations of the OM data when 
plotted in Fig. 5. In approximately 50% of the sources, (typically those for which x 2 < 13) the fits have a probability > 0.05 that the 
model is an acceptable representaion of the data. The second line of data for each source gives the fractional mass accretion rate (Kerr 
models), fractional luminosity (Schwarzschild models) and the cosine of the inclination angle - the models sampled the full parameter 
space of these parameters and we include their values for completeness. N.B. In order to determine errors (to within 99% confidence) 
on the resultant masses, we have obtained the volume in parameter-space that is contained within Ax 2 of the minimum x 2 for each fit. 
This does not enable us to constrain the values of mass accretion rate, luminosity or inclination angle but gives us approximate error 
ranges for the masses. These error ranges should, however, be used with caution since they do not reflect exactly those points which were 
included within the A% 2 + Xmin vomme of parameter-space. 
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Figure 6. Comparison of the XMM photon indices with those 
of ROSAT (Laor et al. 1997; the energy indices quoted have 
been converted to photon indices) and ASCA (George et al. 2000; 
Reeves & Turner 2000). The dotted lines represent the point at 
which each pair of telescopes would obtain equivalent results. 



The BQS sample has also been observed previously 
at harder energies using the ASCA satellite. George et 
al. (2000) found a range F = 1.5-3.0 which is consistent 
with the XMM-Newton values for the Laor et al. sample. 
PG 1411+442 was the only source in the George et al. (2000) 
sample that could not be modelled successfully with a single 
power law in this range and, similarly, we have found that it 
warrants more detailed fitting than performed in this paper. 
ASCA data from a larger sample of 62 quasars was studied 
by Reeves & Turner (2000) and similar values obtained for F 
in the 2-10 keV range. Comparison of the XMM and ASCA 
photon indices is shown in the bottom panel of Fig. 6. 

Table 8 provides a summary of X-ray properties of in- 
dividual sources. In only four of the quasars studied in 
this paper - PG 1216+069, PG 1309+355, PG 1425+267 
and PG 1440+356 - could an additional neutral absorp- 
tion component be detected. Warm absorbers have been 
discovered in the XMM-Newton spectra of PG 1114+445 
and PG 1309+355 (Ashton et al. 2004) and possibly 
PG 1425+267 and PG 1411+442 (Ashton, Brocksopp et 
al. in prep.). These findings are consistent with Laor et 
al. (1997), who identified a distinct group of "X-ray-weak 
quasars" consisting of PG 1001+054, PG 1411+442 and 
PG 1425+267 on account of their absorption properties. 
We note that these same X-ray-weak objects are among the 



sources which show signs of UV absorption in their spectra 
(Laor & Brandt 2002) Reeves & Turner (2000) found absorp- 
tion in just over half of their sources, mainly in the more dis- 
tant objects but also in some of the low-redshift radio-loud 
quasars. In particular they discuss possible warm absorber 
components in PG 1226+023 and PG 1425+267, again con- 
sistent with our results. This fraction of Laor sources with 
a warm absorber is notably less than the 50% of a sample 
of Seyfert 1 galaxies (Reynolds 1997), as discussed by Laor 
et al. (1997) and Ashton et al. (2004). 

We find little evidence for iron lines in the XMM- 
Newton spectra although the upper limits to the equiva- 
lent width allow for appreciable line emission in some other 
sources. Detections for which the equivalent width could be 
constrained (out to 99% confidence) were obtained for just 
twelve sources but to low significance (< 2a). Our upper 
limits are consistent with the ASCA results since Reeves & 
Turner (2000) report iron line detection in PG 1114+445, 
PG 1116+215, PG 1226+023 and PG 1425+267. We fur- 
ther note that Page et al. (2004a) detect weak iron emission 
in PG 1226+023 when they combine a series of nine obser- 
vations together. Similarly our results are consistent with 
those of Porquet et al. (2004). 

On fitting a power-law to the 2-5 keV range for each 
source we find that there is a soft excess in each case where 
we do not find evidence for a warm absorber. This is in con- 
trast to George et al. (2000) and Reeves & Turner (2000) 
who find a soft excess in just a few of their low-redshift 
quasars, although this is likely to be a reflection on the cal- 
ibration of ASCA at low energies. Porquet et al. (2004) find 
more evidence for a soft excess component but, as also dis- 
cussed by Page et al. (2004b), the derived black-body tem- 
peratures are thought to be too hot to be direct emission 
from a thin accretion disc. Instead they consider that UV 
disc photons are upscattered to X-ray energies via a Comp- 
tonizing medium. This is consistent with the correlations 
we find in Section 4 - the strong correlation between the 
soft and hard power-law photon indices would suggest that 
the 0.3-10 keV emission is all produced by the same Comp- 
tonizing medium and not directly by the disc. In this case 
the optical/X-ray continuum luminosity correlations would 
also be expected since an increase in accretion rate could 
lead to an increase in luminosity at all frequencies. 

It should therefore be possible to fit a Comptonized 
disc model to the optical-UV-X-ray region. In Section 5 
we attempted to do this using the Czerny & Elvis (1987) 
model (both Kerr and Schwarzschild) . While our fits would 
be much better constrained if we were able to sample the 
peak of the emission (presumably in the EUV range), some 
of the published mass estimates produce curves which are 
considerably less luminous than the observed optical/UV 
data. 

Aside from the mass estimates, however, it is clear that 
the accretion disc models presented here do not account for 
the soft X-ray excesses of these low redshift quasars. This 
result, coupled with the correlation between the soft and 
hard X-ray slopes, leads us to question what sort of phys- 
ical models are able to account for the observed spectra. 
If, as the tight correlation would suggest, the soft and hard 
X-ray parts of the spectra are produced via the same emis- 
sion mechanism and in the same emitting region, it is likely 
that we are looking for a non-thermal mechanism, such as 
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Table 8. Summary of the X-ray properties of each individual PG quasar. 
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those invoking Comptonization and/or synchrotron, as ap- 
pears to be the case for the broad-band spectra of X-ray 
binaries in the low/hard state (e.g. Markoff et al. 2001). 
Models invoking thermal blackbody and/or bremsstrahlung 
components are non-physical at high energies and so the cor- 
relation would suggest that they are equally inappropriate 
at low energies. A two-temperature disc-corona model (such 
as Kawaguchi, Shimura & Mineshige 2001) requires the low- 
and high-energy power-law components to be produced in 
different regions and via different mechanisms, creating dif- 
ficulty in producing such a tight correlation. 

The failure of the accretion disc models to fit the soft 
X-ray spectrum also poses a possible problem for models 
of Narrow Line Seyfert 1 (NLS1) galaxies which invoke an 
accretion disc with large inner radius for the soft X-ray 
excess (e.g. Pounds, Done, Osborne 1995). With H/3 line 
widths of 1000-2000 kms" 1 , PG 1115+407, PG 1440+356 
and PG 1543+489 can be classed as NLS1 objects; we find 
that these sources are more variable and have steeper spec- 
tra than most of the other sources in the sample, consistent 
with other NLS1 objects. However, the models yield consid- 
erably larger masses than the 10 6 -10 7 M© that would be 
expected from NLS1 models. Alternative models for these 
objects might instead invoke a more face-on inclination an- 
gle and/or a more distant broad line region (e.g. Boiler et 
al. 1992, Wandel 1997, Puchnarewicz et al. 1992). 

We should question whether more complicated multi- 
wavelength models are required. It is known that, like many 
quasars and active galaxies, PG 1226+023 (3C 273) has jets 
which emit at radio, optical and X-ray wavelengths (e.g. 
Courvoisier 1998 and references therein) and it is likely that 
the other radio-loud PG quasars have similar structures. 
Some radio-quiet quasars have now had their jets imaged at 
radio frequencies (e.g. Blundell et al. 2003a) and, although 
none of the Laor radio-quiet objects have been resolved, they 
arc emitting at significant flux densities of at least a few mili- 



Janskys. Fitting disc models to any source, for which a jet 
is making significant contributions to the optical or X-ray 
luminosity, will be meaningless. 

Four of the sources studied here have a significant 
number of radio points from which we can determine the 
spectral index; of these two are radio-loud (PG 1425+267, 
PG 1512+370) and two are radio-quiet (PG 1116+215, 
PG 1216+069) 4 . Brocksopp et al. (2004) found that for 
the radio-loud high redshift quasar [HB89] 2000-330 the 
spectral index of the X-rays and radio emission were con- 
sistent and suggested that the X-ray emission might share 
a common origin with the radio, such as synchrotron self- 
Comptonization of jet photons. We see a similar relationship 
in just one of the low redshift radio-loud objects studied in 
this paper, PG 1512+370. In contrast, for the two radio- 
quiet quasars studied here we find that the spectral indices 
of the X-rays are steeper than those of the radio. Conversely, 
the radio-loud quasar, PG 1425+267, has a steeper radio 
spectrum than X-ray. The radio spectral indices were in the 
range —0.6 to —0.9, consistent with optically thin emission. 
However, we also note that any curvature in the spectra may 
complicate comparison of the slopes and so a more thorough 
analysis should include study of the X-ray vs. radio flux cor- 
relation. We will address this in a future work. 

For these same four radio-loud sources we also use the 
formulae of Marscher (1983) to estimate the amount of syn- 
chrotron sclf-Compton (SSC) emission expected at X-ray 
wavelengths, assuming that the radio emission was all emit- 
ted in a jet. We find that less than 1% of the X-ray emission 
could be from such a source and so an even smaller frac- 



4 We note however that the concept of "radio-loud" and "radio- 
quiet" may have been introduced artificially through obtaining 
quasar samples according to their optical properties and choosing 
specific radio observing frequencies (see e.g. Blundell 2003b and 
references therein) 
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tion is to be expected in the radio-quiet sources. Therefore 
we must look to alternative sources of X-ray emission and 
more sophisticated jet models in order to account for the 
observed spectra. For example, Zdziarski (1986) presents 
a model invoking SSC to explain the IR-X-ray spectra of 
radio-quiet AGN; similarly the model of Falcke & Biermann 
(1995) produces a significant level of SSC at X-ray ener- 
gies and may be applicable to radio-quiet objects as well as 
radio-loud. Again, we will be investigating these models in 
a future work. 



7 CONCLUSIONS 

We have analysed XMM-iVewton/EPIC X-ray spectra for 
the Laor et al. sample of PG quasars and found that a model 
invoking a broken power-law and iron line is an acceptable fit 
in the majority of cases. Those sources for which the model 
is a poor fit seem to require some sort of warm absorber 
at low energies and one source appears to need a further 
higher energy component. These discrepant sources will be 
studied in more detail in a future paper. Iron lines were de- 
tected in a number of sources but with significance < 2a. 
We have confirmed correlations between various optical and 
X-ray line/continuum parameters; in particular we find that 
the power-law indices of the low and high energy power-laws 
are well correlated, suggesting that both components have a 
common origin. We have also analysed XMM-iVewton/OM 
photometry and combined this with the X-ray spectra and 
archive radio-UV photometry in order to compile broad- 
band spectral energy distributions. Fitting these spectra to 
a Comptonized disc model show that the soft excess appears 
to be independent of the disc model, adding support to our 
interpretation of the correlation between the soft and hard 
X-ray spectra. 
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